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We have analyzed the vahdity of uncertainty relations between the fluctuations of thermodynam- 
ically conjugated extensive and intensive variables within the field of statistical mechanics. Analysis 
is presented for two particular examples of small lattice systems that are in contact with reservoirs 
of comparable sizes: an Ising paramagnet and an Ising ferromagnet. The numerical results enable 
determination of the range of applicablity of the proposed relations. Due to the fact that the exam- 
ples correspond to systems described by discrete variables, the uncertainty relations are not valid if 
the probability of boundary states is too high. 
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I. 



INTRODUCTION 



Uncertainty relations can be established in any physi- 
cal theory involving random conjugated variables charac- 
terized by probability distribution functions. The present 
paper deals with uncertainty relations in statistical me- 
chanics. 

Statistical mechanics combines equilibrium thermody- 
namics with a probability postulate [l|, Following 
the standard Gibbs ensemble theory, in the microcanon- 
ical ensemble it is assumed that all mechanical variables 
are constrained to a fixed value (or are allowed to 
vary within a fixed small range), while their entropy- 
conjugated variables are free to fluctuate. Other ensem- 
bles (Canonical, Grand-Canonical, ...) are introduced by 
relaxing the constraint on one or more mechanical vari- 
ables. This is done by assuming that these mechanical 
variables can be exchanged with a reservoir, which is de- 
fined as a macroscopic body of dimensions much larger 
than those of the system. The exchange is supposed to 
be controlled by parameters (which are characteristic of 
the reservoir) that are identified as intensive variables. 
Formally, the change of variables when passing from the 
microcanonical to other ensembles is undertaken via Leg- 
endre transforms. In the spirit of equilibrium thermody- 
namics, the intensive variables are assumed to have the 
same value both in the reservoir and in the system and 
are usually treated as non-fluctuating variables. There- 
fore, within this basic framework, uncertainty relations 
are not expected to apply [sj. 

While it seems reasonable that intensive variables of 
a very large reservoir (a system characterized by large 
amounts of mechanical variables [l6i |) are treated as non- 
fluctuating variables, there is no reason for such an as- 
sumption when dealing with intensive variables corre- 
sponding to a small reservoir or system. Therefore, we 
will adopt the point of view that, in general, for a system 
in contact with a reservoir (regardless of size) fluctuations 
of intensive variables should occur in the system induced 
by the fluctuations of the corresponding conjugated me- 
chanical variables. From this point of view, uncertainty 



relations are expected to play an important role for a cor- 
rect understanding of the behaviour of equilibrium fluc- 
tuations of thermodynamic variables, with consequences 
in the experimental measurements of these variables. Ac- 
tually this simple picture has received some experimental 
support [3] and thus it is reasonable to ask whether or not 
uncertainty relations can be formulated and what their 
applicability limits are. 



The problem was already discussed by Bohr and 
Heisenberg as soon as Heisenberg proposed the uncer- 
tainty principle in quantum mechanics Q. Therefore, 
during the last century a number of papers were pub- 
lished dealing with this problem. An interesting article 
where these issues have been reviewed is the paper by 
Uffink and van Lith More recently, Velazquez and 
Curilef |2|] revisited the problem within the framework 
discussed some years before by Mandelbrot 0, Gilmore 
and others. The vast majority of these works have 
been formulated from a purely theoretical point of view, 
and applications to particular examples (analytical or nu- 
merical) are scarce. The aim of the present paper is to 
numerically analyze the validity of uncertainty relations 
in two simple lattice models and shed some light on con- 
troversial interpretations [3, 



The paper is organized as follows. In section |TT] we 
summarize a derivation of the uncertainty relations in 
statistical mechanics. In section Hill we numerically ana- 
lyze two examples corresponding to an Ising paramagnet 
and an Ising ferromagnet. Finally, in Section llVI we dis- 
cuss the relevance of the results in relation to tempera- 
ture measurements in small thermodynamic systems and 
draw some conclusions. 
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II. UNCERTAINTY RELATIONS IN 
STATISTICAL MECHANICS 

A. Fluctuations in continuous probability 
distributions 

Following the mathematical treatment proposed in 
Ref. 0, let us consider a continuous probability distribu- 



tion, p (x) dx, defined within the domain f2 £ K", which 
is assumed to vanish at the boundary 9f2. The average 
value of the derivative of any analytical function / (x) 
with respect to one of the components of x, can be ex- 
pressed as: 



dl 
dx, 



g/(x) 

dxi 



= /d"^(^(p(x)/(x))-/(x) 
n 

1 dp (x 



dp (x) 



d"xp(x)/(x) 



p (x) dx., 



dxi 

(/(x)r,(x)) 



(1) 



where we have defined the affinities: 

1 dp{x) dlnp{x) 



^^ (X) 



J3(x) dXi 



dxt 



(2) 



Two interesting particular cases of the result in ([T]) can 
be obtained by choosing / (x) = 1 and / (x) = Xj . In 
the first case one gets 



(r,) = Vi, 

and in the second case 

Let us now define the differences: 

Axj = Xj — (xj) , 



Ar,: = r,. 



(3) 



(4) 



(5) 



(6) 



The variances ((AA)^) and the covariances (AAAB) of 
arbitrary differences AA and AB can be interpreted as 
a scalar products in a vectorial space where the metric 
is given by the probability distribution. In particular, 



the square root of the variance y ((AA) ) of a random 
variable A (which we will denote as the fluctuation of A) 
is the norm of the vector AA in this metric. Therefore, 
in our case, taking into account the Schwartz inequality, 
we can write: 



:Ax,f )((Ar,on>(Ax,Ar,) . (7) 



From ^ and the right-hand-side term can be ex- 
pressed as: 



{Ax, AT,) = {xjT,) - {x.i) (F.) = 4 



(8) 



Thus, we obtain an inequality relating the variances of 
the variable Xi and the affinity Tj : 



{Ax,f) ((AF,;)') > 4 



(9) 



It is worth noting that other authors have used different 
mathematical approaches [^, [ll| to treat this problem. 



B. Application to Statistical Mechanics: 

In this subsection we will discuss how to apply the 
results obtained above in order to establish uncertainty 
relations in statistical mechanics. Actually, the key point 
will be to find a convenient definition of the set $7, the 
probabilities p(x) and the affinities F^. 

We proceed by dividing an isolated universe into two 
interacting parts that can exchange one or more mechan- 
ical variables. One part (usually the smallest one, which 
we wish to study) will be called the " system" , while the 
other part will be the " reservoir" . For the closed universe 
we assume Boltzmann equiprobability and therefore its 
entropy will be given by: 



S/k = \xiW , 



(10) 



where W is the total number of microstates of the uni- 
verse and k is the Boltzmann constant. 

Let us consider the set of mechanical variables de- 
scribing the system X = (Xi, X2, • • • X„). The cor- 
responding variables for the universe X'^ are constant 
(closed). Therefore, for the reservoir we have X^ = 
X*^ — X. From knowledge of the number of microstates 
of the universe VF(X) for which the system variables take 
values within (X, X -|- dX) we can determine the proba- 
bility density: 



p (X) dX 



W (X) dX 



(11) 
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where Qq is an elementary volume in that will be ir- 
relevant for further developments. It is introduced here 
in order to have consistent units. Assuming that this 
probability density p(X) is a finite function, Eq. (jlip can 
be written in terms of a universe dimensionless entropy 
which depends on the system variables: 



and 



(X) = In [W (X)] 



Therefore we have 



p (X) dX = 



(12) 



(13) 



which is none other than Einstein's formula for fluctua- 
tions [13: [3 ■ this point, without making any assump- 
tions about the separability of the function (X) , we 
can define the affinities between the system and reservoir 
as the partial derivative: 



dS^ (X) 
dX, 



(14) 



Note that Ti has units inverse of those of Xi since 
is dimensionless. In some cases, it can be assumed that 
(X) can be split into a part related to the entropy of 
the studied system, S, and a part related to the reservoir, 
S^, so that 



S" (X) = S* (X) + (X^ - X) 



(15) 



Since functions S and only depend on the internal 
structure of the corresponding parts, it then seems natu- 
ral to define intensive variables for the system and reser- 
voir in analogy with thermodynamics: 



dX, 



(16) 



Therefore, assuming separability, the affinity can be ex- 
pressed as: 



It 



(17) 



The results derived in subsection III Al can now be writ- 
ten in terms of the preceding definitions. Without assum- 
ing separability, from Eq. ^ we get: 

(r«) = o, (18) 

and from the inequality relation in ([S]), we can write: 



y((AA,f)(rf)>4 



(19) 

When separability holds, equations (fT8)) and (fT9)) become 
(7f-7.)=0 (20) 



((AX, 



[if 



It 



> 5,. 



(21) 



respectively. These are the main results that we will nu- 
merically test in the examples in section IIIII Note that 
for the particular case of choosing energy as a mechanical 
variable, in the separable case, the above two equations 
read: 



and 



- /?) = 



(Ai?2)((/3fl-/3)2^ >1 



(22) 



(23) 



where we have used the more common symbols /? and 
to represent the intensive variables conjugated for the 
energy of the system and reservoir. 

Equations ([T51) . (pn| or are essentially probabilis- 
tic versions of the Zero-th Principle of Thermodynam- 
ics [12]. Because of the fiuctuations of the system and 
reservoir, the equilibrium conditions are not an abso- 
lute equality, but only apply to statistical averages of 
the intensive variables. It is important to note that with 
the separability assumption the conjugated variables of a 
system appear as functions of extensive variables of the 
same system (as in thermodynamics) and that the same 
happens when the roles of the system and reservoir are 
interchanged. 

Equations (jl9p . (PT|) or state that fluctuations 
in conjugated thermodynamic variables are related in a 
form which is similar to the uncertainty relations in quan- 
tum mechanics. These results are usually referred to 
as uncertainty relations in statistical physics and, even 
though their interpretation in quantum mechanics is dif- 
ferent, they both come from the same mathematical base, 
as noted in Ref. 9. 

At this point it is convenient to point out that the pre- 
vious results were obtained after neglecting the contribu- 
tion from the boundary dVL in Eq. ([T]) and after consider- 
ing a continuous character for the variables Xi describing 
the system. Therefore, the uncertainty relation dT51) and 
the equilibrium condition (|T8l) may not hold if states with 
extreme values of the variables Xi (at the boundary dVt) 
are likely to occur. Such states will be called boundary 
states and in the following examples, that have a discrete 
character, it will be important to evaluate their proba- 
bility and check under which circumstances they become 
too high. 



C. First-order approximation 

Let us consider a universe in which the dimensionless 
entropy can be separated between the system and 
reservoir and we focus only on the energy of the sys- 
tem £' as a mechanical variable. Let us assume that 
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or equivalently, the continuous probability density 
p{E)dE, is centred about a maximum at E'^'^. We also 
require that the function has no other maxima, which 
is a typical situation in statistical mechanics, far from 
phase transition points. Exactly at the maximum we 
have r = 13^ — /3 = (otherwise it would not be ex- 
tremal). Fluctuations of F about zero are related to how 
much the system fluctuates near the point i?*^'. 

Fig. [1] shows a schematic representation of this situa- 
tion. The plot above shows the probability of the energy 
fluctuations. The entropies of the system, the reservoir 
and the universe are shown in the middle diagram. At 
the bottom we show the derivatives of the entropies corre- 
sponding to /?, [3^ and F = — /3. The arrows indicate 
the width of the fluctuations of the energy and the in- 
verse temperatures. If the probability density becomes 
narrower, energy fluctuations will decrease, but will 
have a more pronounced maximum. This means that 
the slope of f3^ — /3 (blue line in the bottom diagram) 
increases and consequently the affinity will have larger 
fluctuations. If energy fluctuations are small enough we 



first-order approximation gives: 



E) 



dE 

dp{E) 



{E"" - E) ■ 



(24) 



dE 



{E - E""^) = ^AE , 



where we have taken into account the fact that E^ is 
fixed and we have defined: 



dE^ 



dp'^ [E) 



dE 



d(3 {E) 



dE 



(25) 



E<'i 



which plays the role of an inverse joint heat capacity. 
Indeed the plots shown in Fig. [1] correspond to such a 
first-order approximation, where the probability is Gaus- 
sian, the entropy densities are parabolas and the inverse 
temperatures and affinity behave linearly. 

Note that within this first-order approximation, con- 
dition (HU yields {E) = E'"i. Besides, the fact that 
(3^ — P (X AE makes (^5)) an equalty (it is indeed a prop- 
erty of the Schwartz inequality). Thus, replacing ((24)) 
into we find 



s 




CO. 
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p^-p 











System energy E 



FIG. 1: Schematic representation of the dependence on the 
system energy E of the probability density (above), the en- 
tropies 5''^, and S (middle), and the inverse temperatures 
P, and the affinity F = - /3 (below). The arrows in- 
dicate the width of the typical fluctuations of the difi'erent 
variables. 

can expand the function r{E) about the point E'^'^. A 



Hence 



{AEr) = 



{AEY 



dl3^ [E) 



I. 



dE 



dp [E) 



dE 



(26) 



(27) 



If we define the system heat capacity C and the reser- 
voir heat capacity as 



dp_ 
'dE 



1 



1 dT 



kT^ dE kT^C 



dp' 



dER 



dp' 



dE 



1 



1 dT 



kT^ dE kT^C^ ' 



(28) 



(29) 



where T ^ {k (/?)) ^ stands for the equilibrium temper- 
ature, we get: 



(^{AEf 



kT'^C'^C 



(30) 



which is essentially the same result obtained by Lindhart, 
as discussed in Ref. 6;, but seen from within a different 
framework. At the same order of approximation we can 
compute: 



(F^) 



1 



{AEY 



kTKJC^ 



c 



kT^C^iC^ + C) ' 



(31) 



(32) 
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ami 



kT'^C{C^ + C) ' 



(33) 



((A/3A/3«)) = V(A/3«)^ {Apf = ^^^^ 



1 



(34) 

From the previous equations, for an infinite reservoir 
(C^ ^ C) we obtain the expected results 



(35) 
(36) 
(37) 
(38) 
(39) 



((A/3^)')=0 



((A/3)') = ^ 



fcT^C 
((A/3A/3«))=0. 



This definition is not separable as a product of terms for 
the system and the reservoir, not even when the probabil- 
ity factorizes. As a consequence, the concept of intensive 
variables cannot be thought of as a property of a single 
system. Instead, in this case, we shall deal with affinities 
that are a property of each pair of systems in equilib- 
rium. In the first example in the next Section Hill we will 
compare the results obtained with the two prescriptions 
in order to evaluate how much the choice can affect the 
results for small systems. 



III. 



NUMERICAL SIMULATIONS 



In order to discuss the applicability of the above ex- 
pressions, we will test them in two numerically treat- 
able examples: a non-interacting lattice system (Para- 
magnetic Ising system) and an interacting lattice system 
(Ferromagnetic Ising system). We will focus on the valid- 
ity of expressions (jlSp and (jl9p for various universe sizes, 
various system sizes within each universe and different 
values for the extensive fixed variables of the universe. 
Our aim is to deepen our understanding of the validity 
of these equations. 



D. Applicability to discrete systems 

Even though many macroscopic systems can be stud- 
ied as continuous, their nature is usually discrete, and, 
therefore, equations and ((TO)) should be seen as ap- 
proximations. We will now discuss how to deal with dis- 
crete systems. 

The main difficulty in the application of equation (jl9p 
to discrete systems, is how derivatives should be taken in 
the definition of intensive variables. We must consider, 
not only the problem of finding whether defining forward, 
backward or central derivatives can affect the result, but 
we must also take into account the fact that 4^h\p and 
P~^lxP ^^^'^ different corresponding discrete expressions. 
Let us discuss two possible definitions. 

If we consider the second part of Eq. ([5]) a possible 
(centred) definition of F is 



Inpi+i - lnp,_i 



(40) 



which is separable if the probability p factorizes. Never- 
theless this prescription also raises some problems: since 
the logarithm goes to minus infinity as probability goes 
to zero, one must be really careful at the edges of the do- 
main and change the centred prescription to backwards 
and forward prescriptions. 

An alternative definition of F^ is the one suggested 
from the first part of Eq. ([2]): 



F,; 



Pi+1 -Pi-l 



(41) 



A. Paramagnetic Ising system 

Let us firstly analyze the validity of the formalism 
developed above for a model without interaction. We 
shall consider a universe consisting of non-interacting 
spins Si that can take two values ±1, depending on 
whether they are parallel or antiparallel to an external 
field H (which is a fixed parameter of the model) . The 
Hamiltonian is given hy % = — X]i=i SiH . Since the uni- 
verse is isolated (constant energy) the number of spins 
in the -1-1 state and the number of spins N'^ in the 
-1 state (fulfilling N^l + N^ = N^) are fixed. Therefore, 
the number of possible configurations of the universe is 
given by: 



W 



(42) 



We now consider a partition of the whole universe into 
the two sets: the system and the reservoir so that — 
N^+N. We will choose the number of spins in the +1 
state in the system as the internal variable of interest in 
this example (this choice gives simpler expressions than 
choosing the energy). Its probability function will be 
given by: 



PiN+) 



m 



N+\{N ■ 
N)l 



N+y. 



{N^~N+)\{NU 



N -N^ 



(43) 
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where , and N are parameters. Note that the 
variable ranges between the extremes corresponding 
to its minimum value given by iV""" = max(0,-/V^ — 
+ N) and its maximum value given by iV™°^ — 
mm {N^,N). Using the probability function (gS]) we can 
explicitly compute the average value (A^+) and its vari- 
ance as 

((AA^+)^) = {(N+f) {N+r (44) 
For illustrative purposes, in this first example we com- 



N^=18 



N^=28 



' N=l 

N=2- 
N=3 - 
N=4 
N=5 
0.8 % N=6 
|\ N=7 
^ l\ N=8 
^ " N=9 



+ 



2^0-4 



0.2 



i ■ ? 




'0 3 6 9 12 15 4 8 12 16 20 24 28 




N 



12 15 



8 12 16 20 24 28 



FIG. 3: Test of the validity of the equality of the (7) and 
{7^) as a function of for a paramagnetic system, using 
the prescription from Eq. (|40l) 



first proposed option, according to Eq. PO)) . enables us 
to benefit from the separation of the entropy density. As 
the spins are non-interacting, when we define the entropy 
using the logarithm and neglect irrelevant constants, it 
can be trivially separated as: 



FIG. 2: Probability of the boundary states p(iV^'") + 
p{N^°'^) as a function of A'^+ for a paramagnetic system for 
different universe sizes and system sizes as indicated by the 
title and legend respectively. 

pute the probability for the system to reach boundary 
states p(Af™™) + piN"^"'-"). Fig. [2] shows the value of 
this probability for two universes of sizes — 18 and 
= 28 as a function of and for different values of 
the system size N. This result shows that the probabil- 
ity of the boundary states departs from zero when the 
system or reservoir are too small and when the value of 
iV^ is too close to its extremes and this fact constrains 
the number of available states in the system. 

So far there has been no need to take one of the pre- 
scriptions described in section III Dl However, in order 
to perform the study of the conjugated intensive vari- 
ables we must now choose one of the prescriptions. The 



S{N+) = lnj3(iV+) = SiN+) + S'^iN^ - N+) (45) 
with 

and 

_ (N'^ -Ny. (47) 

^ ^ {N^ - N+y.{NU - N - N'l + N+y. 

To estimate the intensive variables of the system (7) and 
reservoir (7^) we approximate the derivatives of S and 

by a finite difference scheme. Since the variable 
is bounded, we will use central differences everywhere 
except at the extremes, where we will use forward or 
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FIG. 4: Test of the validity of the uncertainty relation as 
a function of A''+ for a paramagnetic system, using the pre- 
scription from Eq. (|4Up . The red line shows the proposed 
lower bound. 



backwards differences. Thus: 




{N - N+ + 1){N - N+) 



(48) 



{N+ + 1)N^ 



and 



S^iN+ + 1) - S^{N+ - 1) 



In, 



- iV - 7V|^ + iV+ + \){NV -N - NY + N+) 



(iVy-iV+ + l)(7V^-7V+) 



(49) 



Fig. [3] shows the difference (7-'' — 7) as a function of N^, 
for universes of sizes = 18 (left panel) and = 28 
(right panel) and different system sizes N (ranging from 
1 to /2). Notice that the roles of the system and 
reservoir are, in this case, symmetric. As can be seen, the 
difference (7^ — 7) vanishes when iV^ (which determines 
the available energy) is not close to its upper or lower 
extremes. It is clearly observed that the range of validity 
of the result in Eq. (I^IH) increases when the system size 



FIG. 5: Test of vanishing of the average affinity (F) for a 
paramagnetic system, using the prescription in Eq. H4ip 



N approaches half the universe size since we have 

more available states and therefore the importance of the 
boundary states on the computed averages decreases. 

Fig. U shows the test of the uncertainty relation (pij) 
as a function of the total number of up spins iV^^. For 
= 18 and = 28, and different values of N. Sim- 
ilar conclusions to those observed in Fig. [3] can be de- 
duced: when N and iV^ are far from their boundary val- 
ues, the uncertainty relation is fulfilled. Only when the 
fluctuating variable is strongly constrained to a tiny 
number of states by the system size N or the available 
fixed universe energy 7V^, do the surface terms neglected 
in the derivation of the uncertainty relation play a role 
and the bound is broken. 

Alternatively, we can study the second method pro- 
posed for discrete systems and define the joint affinity F 
as in Eq. (j4ip . In this case, however, it is not possible 
to separate F into two terms, one corresponding to the 
system and another one corresponding to the reservoir. 
Numerical computation for the same cases studied be- 
fore allows us to plot the results shown in Figs. [5] and 
[51 The same qualitative results concerning the validity 
of and can be concluded. Note that the values 
of (F) in Fig. [5] are smaller than those in Fig. [31 There- 
fore, for discrete systems, the choice of the approximation 
method discussed in section III Dl may affect the results 
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FIG. 6: Test of the validity of the uncertainty relation for a 
paramagnetic system using Eq. H41[) . The red hne shows the 
proposed lower bound. The symbols correspond to the values 
of A'^ indicated in the legends of the previous figures. 



FIG. 7: (F) plotted against the energy of the universe for a 
15x15 and a 8x8 universe, and for various sizes of the corre- 
sponding systems. As the size of the system increases, the 
expected value of F goes to zero, as expected. 



from a quantitative point of view, even in such a simple 
case. 



B. Ferromagnetic Ising system 

In this second example we study the case of a discrete 
system with interactions. Since interaction can induce 
correlations, in this case we should not assume that the 
entropy density is separable and thus we should re- 
strict ourselves to only testing the validity for equations 
([T5)) and . Let us consider a universe consisting of an 
Ising model on a two- dimensional square lattice with size 
— ly X LF ^ with nearest- neighbour interactions and 
periodic boundary conditions. The system of interest will 
be a rectangular subset of spins with size N = x Ly. 
Even though a general treatment requires the simultane- 
ous analysis of energy and magnetization fluctuations, in 
the absence of an external field, it is enough to analyse 
only the energy. The Hamiltonian of the universe reads 

n.n. 

= - ^ s,sj , (50) 



where Si = ±1 are the spin variables and the sum extends 
over all the 2N^ nearest-neighbour pairs of the universe. 
Note that the energy of the universe will be discretized 
in steps SE^ = 4 and ranges from —2N^ to 2N^ For 
any given configuration, the evaluation of the energy of 
the system will be done by considering that the energy 
associated with the bonds along the boundary of the sys- 
tem is equally divided between the system and reservoir. 
Therefore, the energy of the system is discretized in steps 
SE = 1. 

In order to perform a microcanonical analysis of the 
averages at constant E^ we have used the following pro- 
cedure. We have randomly generated spin configurations 
and classified them according to the universe energy. In- 
side the set corresponding to the same value of E^ we 
have computed the fraction of configurations that cor- 
respond to every possible value of the energy E of the 
system and performed a statistical estimation of its prob- 
abihty p{E). 

It is important to note that, because of the finite size 
of the random sample of states that we have generated 
(3 X 10^ configurations), we do not obtain states with an 
a priori small but non-zero probability. This means that 
there is a certain number of possible values of E for which 
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N^=8x8 



n"=15x15 



1.5 h 



0.5 
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FIG. 8: Test for the fluctuation theorem in discrete systems. 
Symbols are the same as in the previous figure. As in figure 
[21 the red line shows the proposed lower bound. 



we wrongly estimate zero probability. As the logarithmic 
prescription (|40| discussed in IIIDI would give problems 
with infinities at the neighbouring values, we have only 
used the second prescription given by Eq. (1411) . Even 
in this second case, in order to avoid problems with the 
computation of the variance of F (for which a pi survives 
in the denominator), the states with an estimated zero 
probability have been excluded: it is as if they could not 
be reached by the system. 

Figures [7] and [5] show the results obtained for simu- 
lations of a 8 X 8 and a 15 x 15 universe, for various 
sizes of the system. These are qualitatively compatible 
with those obtained for the system without interactions 
in Sec. IIII Al The uncertainty relations are fulfilled ex- 
cept in the very extreme case in which the number of 
available states of the system and/or reservoir is so small 
that the boundary effects are dominant. 



IV. DISCUSSION AND CONCLUSIONS 

We have shown that uncertainty relations can be es- 
tablished in statistical mechanics when considering a uni- 
verse (isolated from the outer world) separated into two 
parts, called the system and reservoir, which are in ther- 



modynamic equilibrium and interact by exchanging "me- 
chanical properties" such as energy, volume, particles, ... 
If Xi is one of these properties corresponding to the sys- 
tem, we have found that the following inequality must 
be satisfied, y/{{AXi)^)({ri)^) > 1. Here, F, is an affin- 
ity conjugated to the variable Xi which determines the 
departure from strict equilibrium of actual microscopic 
configurations of the two interacting systems and which 
satisfies (F^) = 0. Furthermore, when the entropy sat- 
isfies separability, these affinities simply reduce to the 
difference 7/^ — ji of the intensive variable conjugated to 
Xi, in the reservoir and the system. Assuming that the 
intensive variable of the reservoir does not fluctuate as 
expected for very large reservoirs, the proposed inequal- 
ity reduces to the more standard uncertainty relation in- 
volving only system variables 



V((AX,)2)((7,)2) > 1 



(51) 



Furthermore, assuming that the fluctuations in the sys- 
tem are small, by expanding to first order, the uncer- 
tainty relation leads to the following equality: 



{{AX, 



1 



(52) 



It is interesting to point out, than within this ffist-order 
approximation, when the relevant conjugated variables 
are the energy and the inverse temperature, combining 
Eqs. dSni) and one gets: 



V((Ai?)2)((A/3)2) = C^/{C + C") < 1 



(53) 



where C and are the heat capacities of the system 
and reservoir. Therefore, in this last case fluctuations of 
inverse temperature are not, in general, expected to be 
proportional to energy fluctuations [4] . Only in the limit 
C/C^ < 1 is the standard result ^ {{AEf ) {{A^f) = 1 
recovered. 

The derivation of the above results has required the 
assumption that the probability for the boundary states 
(on the boundary of the set of mechanical variables Vl) 
vanishes. This assumption, which is reasonable for uni- 
verses described by continuum variables, might be incor- 
rect for the discrete case, when universes are small or too 
constrained. 

We have numerically studied the validity of the above 
results in two different mechano-statistical lattice sys- 
tems: the Ising paramagnet and the Ising ferromagnet. 
We have found that the uncertainty relations hold ex- 
cept in the case in which the boundary states have a high 
probability of occurring. In particular, from the exam- 
ples in section Hill we can numerically evaluate a bound 
for this probability. Fig. IH] shows the uncertainty prod- 
uct ((AXi)2}(r|) against the probability of the bound- 
ary states for the Ising paramagnet (with the two possi- 
ble definitions of the affinities for discrete systems given 
by Eqs. (HUl) and (jH])) and the Ising ferromagnet. The 
figure clearly shows that the breakdown of the uncertain- 
ity relation may occur only when the probability of the 
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boundary states is above 0.01. In all the other cases, the 
uncertainty product is larger than one and, in fact, tends 
to 1 from above when the size of the universe increases, 
as can be seen in the inset of FiglHl 

These results are relevant from an experimental point 
of view. Focusing on the case of the energy E and the 
inverse of temperature /3 as conjugated variables, we can 
identify the reservoir with a given body and the system 
with the thermometer used to measure its temperature. 
Commonly, the thermometer is chosen so that it has a 
heat capacity much smaller than the heat capacity of the 
body. However, for small bodies, when the heat capacity 
of the thermometer and the body are comparable, the 
results derived here must be taken into account: (i) the 
inverse temperature (/3) measured by the thermometer 
con strongly fluctuate, but still its average value must 
be identified with the average of the inverse temperature 
of the body {(3^); and (ii) the fluctuations of the differ- 
ence /3^^ — (3 should be related to the energy fluctuations 
through the proposed uncertainty relation. 

Future extensions of this work will deal with the anal- 
ysis of systems described by continuous variables, where 
direct counting of states is not straightforward. The re- 
sults will be presented elsewhere. 



FIG. 9: Uncertainty product against the probability of bound- 
ary states for the paramagnetic Ising system and the ferro- 
magnetic Ising model. Data includes different universe sizes, 
system sizes, and different values of the mechanical variables 
of the universe. We have also included the data correspond- 
ing to the two different methods proposed for the evaluation 
of the affinities for discrete systems, given by Eqs. (140 1) and 
(|41|l . indicated in the legend by the numbers between paren- 
thesis. The inset shows a magnification of the region above 1 
for probabilities below 2 x 10~*. 
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